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Environmental context strongly affects many host-pathogen interactions, but the underlying causes of these 
effects at the individual level are usually poorly understood. The amphibian chytrid fungus has caused 
amphibian population declines and extinctions in many parts of the world. Many amphibian species that have 
declined or have been extirpated by the pathogen in some environments coexist with it in others. Here we show 
that in three species of rainforest frogs in nature, individuals' probability of infection by the amphibian chytrid 
fungus was strongly related to their thermal history. Individuals' probability of infection declined rapidly as 
they spent more time above the pathogen's upper optimum temperature. This relationship can explain 
population-level patterns of prevalence in nature, and suggests that natural or artificial selection for higher 
thermal preferences could reduce susceptibility to this pathogen. Similar individual-level insights could 
improve our understanding of environmental context-dependence in other diseases. 

Many symbiotic relationships can be modified extensively by environmental conditions 1 . Harmless sym- 
biotes can become pathogenic, and mutualistic relationships can break down 1 . The emerging infectious 
disease chytridiomycosis, caused by the amphibian chytrid fungus Batrachochytrium dendrobatidis 
{Bd), is an excellent system for examining such environmental context dependence. Chytridiomycosis has caused 
amphibian population declines and extinctions in many parts of the world 1,2 , yet many amphibian species that 
have declined or have been extirpated by the disease in some environments coexist with Bd in others 1 " 4 . 
Additionally, some infected amphibians can survive for years in the wild with no clinical signs of disease 5 , and 
some species that coexist with endemic chytridiomycosis in the field can be highly susceptible to Bd in the 
laboratory, experiencing 100% mortality 6 . 

It is known that environmental temperature can determine the progress and outcome of Bd infections in the 
laboratory. In culture, Bd grows best between 17 and 25°C 7 . In the laboratory, low and/or fluctuating tempera- 
tures can retard the pathogen's growth, and elevated body temperatures can clear frogs of Bd infection 6 . In nature, 
population-level infection prevalence and host mortality rates are often correlated with ambient environmental 
conditions; the infection rates are highest during cooler months 2 4 and at higher elevations 4 . Similar context- 
dependence occurs in other host-pathogen systems 1 . 

One study 8 has shown that after the first appearance of chytridiomycosis in a susceptible population, mean 
body temperatures of frogs increased, which should reduce the negative effects of the pathogen. Increases in mean 
body temperature at the population level could reflect adaptive responses of individuals to infection, shifting 
thermal preferences to produce "behavioural fever" 8 . However it could also result from selective sweeps in which 
individuals that attain higher temperatures for other reasons 9 are more likely to survive during outbreaks of 
chytridiomycosis. Understanding such population-level responses and determining their implications for the 
management of amphibian populations requires information on the individual-level patterns that underlie the 
population-level relationships. 

In order to examine the relationship between Bd infection status and individual thermal behavior, we measured 
the Bd infection status and body temperatures in nature of stream-associated rainforest frogs of three species that 
have declined to different degrees due to chytridiomycosis in rainforests of northern Queensland, Australia; 
Litoria lesueuri (least affected by declines), Litoria serrata (intermediate) and Litoria nannotis (most affected). 
Frogs were tracked in the summer/wet and winter/dry seasons at a total of four rainforest sites and their individual 
body temperatures were recorded nocturnally and diurnally every 24 h. 

Results 

Infection prevalence did not differ significantly among sites within seasons for any of the three species (Fisher's 
exact tests, all P > 0.05). Infection prevalence appeared to differ among species and seasons (Fig. la). Our 
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Figure 1 | Individual probability of infection by the amphibian chytrid 
skin fungus and individual thermal history in nature, (a) Prevalence of Bd 
infection in tracked frogs by species and gender and by season. Bars 
represent binomial (Clopper-Pearson) 95% confidence limits, 
(b) Boxplots illustrating distribution of percent of frog body temperatures 
above 25 °C by species and gender and by season. Heavy bars represent 
medians, boxes represent 25% and 75% quartile boundaries, whiskers 
represent ranges, and stars are points outside 1.5X the interquartile range 
from the nearest quartile boundary. Comparing (a) and (b) suggests an 
inverse relationship between prevalence and percent of body temperatures 
above 25°C; this relationship was explored statistically in the models 
presented in Table 2. 

ANOVA showed that the proportions of individuals' body tempera- 
tures that were above 25°C (this temperature was chosen as the cut- 
off because in vitro the growth rate of Bd declines sharply above that 
temperature 7 ) was significantly affected by species, gender, season, 
and their interactions (Table 1; Fig. lb). 

Model fitting and selection of models relating species, gender, 
season, and proportion of individual body temperatures above 
25°C produced three candidate models with delta AICc less than 3 
(Table 2). These were averaged to produce a final model that 
included the effects of species, gender, and proportion of body tem- 
peratures above 25°C (Table 2; Fig. 2). The averaged model shows 
that individual probabilities of infection by Bd decreased strongly 
with increasing percentage of body temperatures above 25°C in all 



Table 2 | Logistic regression models evaluated as possible fits to 
infection status (0 = infected, 1 = uninfected) of tracked frogs. All 
possible one, two, and three-variable main-effects only models 
involving species identity (Species), gender of the individual (Sex), 
season (summer/wet season or winter/dry season, Season), and 
the proportion of body temperatures measured for each individual 
that were greater than 25°C (PA25) were considered. Models 
including more than three main effects or interactions were not 
included to preclude overfitting. The three best fitting models, with 
a total Akaike weight of 0.705, were averaged to obtain the final 
model. Predictions of this model appear in Figure 2 
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'coefficients structurally set to zero; f one-tailed P value due to preexisting hypothesis regarding 
effect direction (higher PA25 should decrease P(infection)). 



three species (Fig. 2). Males and females of the least vulnerable spe- 
cies (L. lesueuri) had the lowest probabilities of infection. Females of 
the two more vulnerable species had moderate probabilities of infec- 
tion, particularly if they rarely elevated their body temperatures 
above 25°C, and males of those species had the highest probabilities 
of infection. Frogs of all species and genders were very unlikely to 
carry Bd infections if 75% or more of their body temperatures were 
above 25°C. No separate effects of season appeared in our best fitting 
models, suggesting that seasonal differences in infection prevalence 
are explained by effects of season on individual body temperatures. 



Table 1 | Analysis of variance for the effects of species, gender 
(sex) and season (cold, dry or warm, wet, season) on individual 
thermal history, measured as proportion of individual body tem- 
peratures that were above 25°C 


Effect 


SS 


df 


MS 


F P 


Species 
Sex 
Season 
Error 


0.740 
0.125 
3.962 
2.920 


2 
1 
1 

97 


0.370 41.262 <0.001 
0.125 4.461 0.037 
3.962 141.107 <0.001 
0.028 



Discussion 

Our study is the first demonstration that individual thermal histories 
affect the probability of infection by Bd of frogs in nature, and one of 
the first examining how individuals' environmental histories affect 
their vulnerability to disease in wildlife. We found that individual 
probabilities of infection by Bd decreased strongly with increasing 
percentage of body temperatures above 25°C in all three species of 
rainforest stream frog examined. This relationship can explain popu- 
lation-level patterns of prevalence in nature 2 4 . 

The effects of species and gender on individual probabilities of 
infection by Bd probably reflect a combination of differences in 
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Figure 2 | Individual probability of infection by the amphibian chytrid 
skin fungus in nature was strongly related to individual thermal history. 

This logistic model, incorporating species identity, gender, and percent of 
body temperatures above 25°C (Table 2), was produced by averaging the 
three best fitted models, each with delta AICC < 3. Color coding 
corresponds to species-level susceptibility to epidemic chytridiomycosis 
(red most, amber intermediate, green least). Predictions of the model 
correspond well (r 2 = 0.739, 4 d.f, P = 0.028) with population-level 
patterns of prevalence when it is used to predict prevalence by season and 
site, factors not included in the model. 

immune function and differences in transmission rates 210 . Our mea- 
sure of thermal history, the percentage of individual body tempera- 
tures above 25°C, is directly related to the known thermal sensitivity 
of Bd in vitro; it thus seems likely that the effects we found were 
caused by direct effects on rates of pathogen reproduction 2,6,7 , how- 
ever they could also reflect responses of frogs' immune function to 
temperature 11 . 

Our findings are not highly consistent with the hypothesis that 
frogs elevated their body temperatures in response to infection 
(behavioural fever) 8 , since the probability of infection declined dra- 
matically with increasing time spent at higher body temperatures. 
However, we cannot entirely eliminate the possibility that infection 
may cause individuals to regulate their body temperatures at higher 
levels than normal. For example, if thermal preferenda vary among 
individuals, we may have only seen infections in those that normally 
have low upper set points and in whom raising the set point has failed 
to elevate temperatures sufficiently to reduce their probability of 
infection. Alternatively, infection could lead to long-term changes 
in thermoregulatory behaviour, so that a propensity to maintain a 
relatively high body temperature reflects what might be termed the 
ghost of infection past. 

Our results may also have a simpler cause. Individuals that choose 
relatively high body temperatures for other reasons, i.e. to aid in 
growth, digestion 9 , or reproduction 3 may have coincidentally de- 
creased their probability of Bd infection. While Bd remains endemic, 
any tendency to maintain higher body temperatures could be rein- 
forced by natural selection. Alternatively, our results could reflect a 
tendency for infected individuals to maintain relatively low body 
temperatures. This could be coincidental; e.g. lower activity levels 
in infected individuals could have reduced their opportunities for 
thermoregulation. It could also reflect manipulation of hosts by the 
pathogen. These hypotheses regarding cause and effect can only be 
resolved by future experimental work. 

Regardless of the cause, our results show that individual histories 
of body temperature can explain a substantial fraction of the sea- 
sonal, elevational, and interspecific patterns of Bd infection preval- 
ence that have been observed in nature 14 . Variation in thermal 
preferences among individuals should provide opportunities for nat- 
ural selection to reduce vulnerability to chytridiomycosis. Mani- 
pulating habitats to increase the availability of warmer temperatures 



may serve as a mechanism to increase the survival of threatened 
amphibians in nature. The effects of global change are likely to vary 
among species and sites; for example increased temperatures might 
tend to decrease the effects of Bd, but increases in cloud or canopy 
cover could tend to increase them. Our results highlight the im- 
portance of the effects of body temperature variation among and 
within individuals on host-pathogen systems involving ectothermic 
animals. 

Methods 

Study sites. The study was conducted at four streamside sites within the tropical 
rainforests of northern Queensland, Australia: Python Creek ("Tully Gorge", 
145°35'E 17°46'S, 200 m asl) and an unnamed creek ("Lower Tully", 145°41'E 
17' 48'S; 70 m asl) in Tully Falls Forest Reserve, an unnamed creek in Kirrama State 
Forest ("Kirrama", 145°52 E 18°11 S; 200 m asl), and Frenchman Creek, in 
Wooroonooran National Park ("Babinda", 145 55' E 17 20' S 20-100 m asl). 
A marked transect at each stream served as a reference for frog locations. 

Study species. Frogs of three species were tracked: the Stony Creek Frog L. lesueuri 
(31 females, 37 males), which has not experienced population declines 12 , the Green 
Eyed Treefrog L. serrata (10 females, 13 males), which declined during epidemics of 
chytridiomycosis in the late 1980s and early 1990s and subsequently recovered 12 , and 
the Waterfall Frog L. nannotis (29 females, 8 males), which was extirpated from many 
sites during chytridiomycosis epidemics 12 . The taxonomy of L. lesueuri was revised in 
2004 13 . Two newly recognized species, L. jungguy and L. wikoxii, occur in sympatry in 
the study region, hybridise, and cannot reliably be distinguished morphologically 13 . 
Our study populations could have included either species, mixtures of both, and 
hybrid individuals. Population declines have not been observed in either species. We 
therefore continue to refer to the study species as L. lesueuri, while recognizing that 
our samples contain unknown proportions of the two morphologically 
indistinguishable species. All three species are large to medium sized hylids 
(males 5.4-12.5 g, females 6.5-41.3 g). 

Tracking. At initial capture, frogs were weighed and swabbed on the ventral surface 
for diagnostic PCR using a sterile rayon swab (Medical Wire & Equipment Co. (Bath) 
Ltd., Wiltshire, UK). Swab samples were analysed using diagnostic PCR or 
quantitative PCR 14 . 

Tracking devices (radio transmitters, or diode tags for harmonic direction finding, 
depending on size of the frog 16 ) were then fitted in situ and frogs were released at their 
point of capture after less than five minutes of handling. Tracking followed protocols 
we have previously published 15 . 

Surveys lasted 16 days and were conducted in the winter/dry season 
(July-September) and the summer/wet season (February- April) at two sites for each 
species. Litoria genimaculata and L. nannotis were tracked simultaneously at the 
same streams during 2004 and L. lesueuri were tracked during 2005. During surveys, 
the location of each frog was determined once during the day (0900-1800 h) and once 
at night (1900-0400 h). Whenever possible, the temperatures of located frogs and 
their substrates were recorded. 

When we could visually locate individual frogs, and were able to reach within 0.5 m 
of them, temperature of the dorsal body surface was measured using a Raytek ST80 
Pro-Plus Non-contact Thermometer (RAYST80, emissivity set to 0.95). This tech- 
nique obtains body temperature readings within 0.5°C of cloacal temperatures 17 . 
When the infrared thermometer could not be used, if the frog was fitted with a 
temperature-sensitive transmitter, we recorded the pulse interval of the telemetry 
signal by timing 90 pulses with a stopwatch. This was converted to temperature using 
individual calibration curves for each transmitter. 

On average, each frog was located 17 times (minimum 5, maximum 28). We 
excluded data from the night following tag attachment from analyses due to the 
potential short-term behavioural effects of handling 18 . Effects of tags on behavior are 
unlikely to persist after the first night of tag attachment 16 . The weights of tracked 
individuals did not change significantly over the study period (Wilcoxon Signed 
Ranks Test; Z — — 1.361, p = 0.173, n = 70). Research was carried out under Scientific 
Purposes Permits issued by the Queensland Parks and Wildlife Service 
(WISP01715204 and WITK01715604) and approved by the James Cook University 
Animal Care and Ethics Committee (A863). 

Data analyses. On swabs analysed using qPCR, we measured between 6 and 10,527 
zoospore genome equivalents per swab. As —20% of our samples were analysed via 
PCR only, we included only infection prevalence in our statistical analyses. Because 
tracking did not occur at the same sites for all three species, we could not include site 
as a factor in orthogonal models. We initially examined our data for possible 
site-specific effects using Fisher's exact tests to determine whether infection 
prevalence differed significantly among sites within seasons for each species. It did 
not (Fisher's exact tests; all p > 0.05). We thus omitted site as a factor in all further 
analyses. We next examined whether individual thermal history, measured as the 
proportion of body temperatures of each individual that were greater than 25°C, was 
affected by species, gender, and season, using an ANOVA in SPSS version 19. This 
analysis showed that thermal histories were affected by all those factors. Finally, we 
examined how Bd infection status (0 — infected, 1 — uninfected) was affected by 
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species, gender, season, and the thermal histories of individuals (proportion of body 
temperatures above 25 C), using a set of models created using the generalized linear 
models procedure in SPSS version 19 with a binary logistic link function (Table 2). We 
used Akaike's information criterion, corrected for finite sample size (AICc) to 
evaluate model fit, using the criteria of Burnham and Anderson 19 . Our best fitting 
models were averaged to obtain the final model 19 . 
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